Europhysics Letters 



PREPRINT 



Restoring a fluctuation-dissipation theorem in a nonequi- 
librium steady state 

T. Speck and U. Seifert 

77. Institut fur Theoretische Physik, Universitat Stuttgart, 70550 Stuttgart, Germany 
PACS. 05. 40. -a - Fluctuation phenomena, random processes, noise, and Brownian motion. 



Abstract. - In a nonequilibrium steady state, the violation of the fluctuation-dissipation 
theorem (FDT) is connected to breaking detailed balance. For the velocity correlations of a 
driven colloidal particle we calculate an explicit expression of the FDT violation. The equilib- 
rium form of the FDT can be restored by measuring the velocity with respect to the local mean 
velocity. 



Introduction. - The fluctuation-dissipation theorem (FDT) relates the correlation func- 
tion of thermally driven fluctuations in equilibrium with the response of the system to a small 
external perturbation [1]. Loosely speaking following Onsager, the decay of a fluctuation is 
independent of whether it has been created spontaneously due to thermal noise or whether is 
has been induced by a small applied force. According to common wisdom, the FDT breaks 
down in a nonequilibrium system. Indeed, a large body of studies deals with violations of the 
FDT in nonequilibrium, in particular, for aging systems (for reviews, see refs. [2,3]), fluids 
and colloidal suspensions driven by shear flow [4,5], and biophysical systems (for an example, 
see ref. [6]). 

An important aspect of nonequilibrium steady states is that detailed balance is broken, 
which requires permanent dissipation of energy. In the case of driven colloidal systems, the 
heat thus dissipated has been identified as "housekeeping heat" [7-9] . Since both the house- 
keeping heat and the violation of the FDT originate from breaking detailed balance, a deeper 
connection between these two concepts can be expected [10]. Since furthermore the velocity is 
the crucial quantity entering the housekeeping heat, this connection should become apparent 
in the violation of the velocity FDT. 

In this Letter, we first derive an explicit expression for this violation. We then demonstrate 
the restoration of the usual equilibrium form of the velocity FDT if the velocity fluctuations 
are properly measured with respect to the local mean velocity rather than with respect to 
the global mean velocity. Our result thus suggests that the decay of fluctuations around the 
steady state is the same, whether they are spontaneously created or externally induced. 

The system. - As a paradigm, we study one-dimensional systems with periodicity I sub- 
ject to stochastic dynamics, where the periodicity is necessary to eventually reach a nonequilib- 
rium steady state (see fig-HJ. Ovcrdamped diffusive motion is then governed by the Langevin 
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Fig. 1 - A colloidal particle moving in a periodic potential V(x) with periodicity I = 2-zr and angular 
position x. A force / can be applied directly to the particle, driving it into a nonequilibrium steady 
state. 



equation [11] 

x(t) = n[F(x)+f p (t)]+r,(t). (1) 

The total force F(x) — —V'(x) + f can be split into the gradient of a periodic potential 
V{x + 1) = V(x) and a non-conservative force /, where the prime denotes the derivative with 
respect to x. The force / p represents a small external perturbation. The thermal noise r){t) 
has zero mean and correlations {r](t)r](T)) = 2DS(t — r). In equilibrium, bare mobility /i, 
temperature T, and diffusion coefficient D are connected by the Einstein relation D — fiT 
with Boltzmann's constant set to 1 throughout the paper. We will keep the Einstein relation 
despite the fact that the system is driven into a nonequilibrium steady state for / ^ 0, which 
expresses the assumption that the driving does not affect the heat bath. The brackets (■) 
represent the average over the thermal noise. 

An unperturbed (/ p = 0) one-dimensional system in a nonequilibrium steady state is 
distinguished from equilibrium by a nonzero, constant current 

j s = / J LF(x)p s (x)-Dp' s (x) ^0. (2) 

We define the stationary distribution p s (x) = exp[— <f>(x)] introducing the generalized potential 
4>{x). In addition, the notion of a local mean velocity 

v a {x) = j s /p s (x) = fiF(x) + D<t>'{x) (3) 

will become crucial. It is given by the drift velocity )J.F(x) plus the gradient of the potential 
4>{x). The local mean velocity ^ s (a;) can be regarded as a measure of the local violation of 
detailed balance. After averaging 

v s = (v s ) = I dx v s {x) Ps {x) - (i (F) - D[ Ps ] l Q = (i (F) = (x) (4) 
Jo 

it becomes the global mean velocity of the particle in a steady state as expected. The boundary 
term vanishes due to the periodicity of the system. 

Fluctuation-dissipation relations. - A small probing force f p (t) is applied after the system 
has been prepared in a steady state described by the stationary probability p s (x). The linear 
response of any observable A to this perturbation is defined as 



RA{t -T) = 



Sf P (r) 



(t > r), (5) 

/ P =o 
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which depends only on the time difference t — r. Due to causality the response function is 
zero for t < r and in what follows we always take t ^ r. We distinguish averages (•} in 
the equilibrium system (/ = 0), averages (■) in the unperturbed system in a noncquilibrium 
steady state, and averages (•) in the perturbed system (/ p ^ 0). 

The various fluctuation-dissipation relations state that the response function of an ob- 
servable is not independent of its correlations. In particular, in equilibrium the well-known 
FDT 

TRf(t-r)=d T (A(x(t))B(x(r))) eq (6) 

relates the response function of A to the correlations of this observable with a unique observ- 
able B conjugated to f p with respect to energy. However, even in a noncquilibrium steady 
state the less known fluctuation-dissipation relation 

TR A (t - T ) = (A(x(t))D<f/(x{T))) (7) 

correlates the observable A to the gradient of the generalized potential 4>{x) [11]. Similar 
relations have been discussed also for chaotic systems [12]. A third relation involving the 
noise can be derived from the Gaussian weight 

^W~ ex p{-i/ drr? 2 (r)| (8) 

of the noise trajectory r)(r) in the time interval to t ^5 t. Since A(x(t)) can formally be 
expressed as a functional A[rf\ of the noise history and variation with respect to /l*/ p (t) is 
equivalent to a variation of t](t) in eq. l|TJl. we have 

TR A (t - t )=d(¥^=dJ [dr,] jjj^m = -Dj [dry] = \ (A(x(t)) v (r)) 

(9) 

through functional integration by parts [3, 13]. 

Velocity fluctuation-dissipation relations. - Because we arc interested in the velocity, a 
fluctuation-dissipation relation involving x instead of the potential <p(x) in eq. JJJ) is called 
for. To this end we insert the Langevin equation Q with the force replaced by the local mean 
velocity into cq. @. In the resulting expression 

TR A (t -r) = (A{x{t))v s {x{r))) - (A(x(t))x(r)) + (A(x(t))r,(r)) (10) 

we use eq. © for the noise term and obtain 

TR A (t -t) = (A(x(t))x(T)) - (A(x(t))u s (x( T ))) (11) 

holding in a noncquilibrium steady state. 

First, we consider the equilibrium case where the position B = x and the probing force f p 
are conjugate. Setting A = x we obtain from eq. 

TRl%t-T) = {x(t)x(r)) eq (12) 

the usual form of the FDT for velocities. We can interchange the order of time derivative 
and taking the average since the probability distribution implicit in the brackets is time- 
independent. 
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In a noncquilibrium steady state, we consequently use cq. (|llfl instead of the FDT JJJJ, 
leading to 

TR ± (t - t) = (±(t)±(r)> - (x(t)v B {x(r))) . (13) 

A somewhat related expression has been derived by a different approach in ref. [10] . Comparing 
cq. H13(l to the equilibrium form (|12|) we identify (£(t)z/ s (x(Y))) as the violation of the FDT 
and define a normalized violation function 

I(t -r) = {[±{t) - v 8 }[u B (x(r)) - v s }) = (x(t)v B (x(r))) ~ v 2 s (14) 

such that I(oo) = 0. 

Restoring the fluctuation-dissipation theorem. - Our crucial observation is that the equi- 
librium form (|12[) can be restored even in a nonequilibrium steady state when we consider the 
relative velocity 

v(t) = x(t) - u B (x(t)). (15) 
Combining cq. Ijllfl and eq. I|13[) . we obtain 

TR v (t -t)= TR±(t - t) - TR Vs {t -r) = («(*)« (r)> . (16) 

This is the main claim of the present Letter. Physically this relation implies that fluctuations 
of the velocity relative to the local mean velocity in a noncquilibrium steady state behave like 
the corresponding response function. In this sense, these fluctuations cannot be distinguished 
from equilibrium fluctuations. 

Connection to the housekeeping heat. - In a nonequilibrium steady state, the heat per- 
manently dissipated in order to maintain the violation of detailed balance is called the house- 
keeping heat [7, 8] . In the time interval to ^ r ^ t it is given as the functional 

QhkMr)] = ii- 1 f dr x(t)p b (x(t)) (17) 

along a single stochastic trajectory x{t). By comparing this expression with eq. (|13f) . a 
connection between energy dissipation and the violation of the FDT for the velocity becomes 
obvious. Indeed, in ref. [9] wc calculated the mean dissipation rate as 

<J = d t (Q hk ) = m- 1 , (18) 

which is constant in a noncquilibrium steady state. For large time differences we get 

(i(tH(x(r))) |f ' rH °° i {vsf = vl (19) 

Since (v^) ^ {^s} 2 , it follows that the violation of the FDT is bounded by the energy dissipa- 
tion rate, 

(x(t)u s {x(r))) ^ [ao~, (20) 

where the equal sign holds for t = r. This inequality specializes the former derivation of 
a general upper bound for the FDT violation [14] to the present case of a driven colloidal 
particle. 

It is interesting to note that we can equate the left hand side of cq. (|20|l involving two 
different times with the mean dissipation rate a by taking into account the change of the 
generalized potential A(j> = <fi(x(t)) — 4>(x(t)) during the time t — r. With 



v s {x(t)) = p s {x(t)) 

I>s(x(t)) P s (x(t)) 



(21) 
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Fig. 2 - Normalized violation function 7(r) over time difference r obtained from 250,000 simulated 
Langevin trajectories for the periodic model potential V(x) = Vbcos2a; (with Vb = \, D = T = 1). 
For large forces /, the function 7(r) shows a non-monotonic decay. 

Fig. 3 - The effective diffusion constant D e ff over the force / for the same potential and parameters 
as in fig. 01 The solid line shows the analytical solution 127H whereas the closed circles are obtained 
numerically from eq. 1241 . The dotted line shows the susceptibility 1251 . 



we see that 



{x(t)v s {x(r))e^) 



[ACT 



(22) 



holds. Since in a nonequilibrium steady state the total heat flow into the heat reservoir is 
Qtot = Qhk — TAcj) [8,9], the term —TAcj> can be regarded as "excess heat" exchanged with 
the reservoir in addition to the housekeeping heat. Its mean value (A<^>) = vanishes and it 
is bounded, whereas Qhk increases in time on average. 

Perspectives. A direct experimental test of the violated or restored FDT, eqs. I|I3|) 
or HI bp respectively, will certainly be a challenge. We therefore suggest first to test an in- 
tegrated version involving the effective diffusion coefficient. As a transport coefficient, it is 
related to velocity correlations by the Grccn-Kubo formula 



D, 



off 



lim 



(x^t)) - (x(t)Y 
It 



dr ([x(r) - v s }[x(0) - v s }) 



Inserting eq. (|13J) leads to 



D 



off 



T 



df P 



f P =o 



/ dr J(r), 
Jo 



(23) 



(24) 



where the first term on the right hand side is the static susceptibility 

d(x) 



POO 

/ AtR ± {t) 
Jo 



df P 



f P =o 



(25) 



and the second term the integral over the normalized violation function I(t). All three quan- 
tities appearing in eq. I)24|) can be measured independently. 
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For a numerical check, we calculate the mean velocity in eq. I|25[) as 




i-cxpH(/ + / p )/r] 



dz exp {±[U{x) - U{x T z)]/T} 



(26) 



with U(x) = V(x) — (/ + f p )x and the effective diffusion coefficient as 




dxl+(x) / dx II 




(x)I_(x), 



(27) 



where now / p = [15]. The violation function I(t) as shown in fig. El is calculated from 
simulated Langcvin trajectories. In fig. |3 we compare the known expression l|27|) with the 
numerical results obtained from eq. (|24|l with excellent agreement. 

Although we discuss in this Letter the one-dimensional case, following our route it is clear 
how to generalize eq. (jlljl to interacting systems with more than one degree of freedom. 
Furthermore, eq. (|16Jl remains valid even if we drop the restriction of the Einstein relation 
in favor of an "effective" temperature T c ff = D(f)//j,, where the strength of the noise D(f) 
depends on the driving. Whether our approach to restore a FDT by referring the velocity to 
the appropriate local mean velocity can be extended to non-velocity like observables involved 
in violations of FDT's remains to be seen. 



Valuable comments by A. Gambassi and an anonymous referee are gratefully acknowledged. 
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